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This talk is a short review on the theoretical issues and uncertainties in the calculation of partial decay rates in 
inclusive B decays. The main emphasis is on charmless semileptonic decays, the B — » X s ^f photon spectrum, 
and the extraction of \V u b\ using model-independent methods. 
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1. Introduction 

Over-constraining the unitarity triangle using many 
different physical processes is a crucial test of the 
CKM picture of CP violation and the validity of the 
Standard Model. The topic of this talk is the theoret- 
ical framework and uncertainty estimation for partial 
rates in B — > X u l~v decays and the extraction of 
the element V U b- Its magnitude is proportional to the 
length of the side opposite the well-measured angle (3 
of the unitarity triangle. 

Because of an overwhelming background from B — > 
X c l~v decays in a large portion of the phase space, 
experimental cuts are employed which suppress such 
background events. The bulk of this talk deals with 
the theoretical techniques used for calculating the sur- 
viving partial decay rates. We may identify two re- 
gions of phase-space for partial rates according to the 
integration domains in the hadronic variables 



P± = E x T \P 



(1) 



where Ex and Px are the energy and 3-momentum 
of the hadronic final state in the B meson rest 
frame. The available phase space in these variables 
is Af2/P_ <P+<M B - 2Ei <P-<M B , where M 2 
is the invariant mass of the lightest possible hadronic 
final state, i.e. the pion mass squared, and Ei is the 
energy of the charged lepton l~ . For the rest of this 
talk we will neglect the finite pion mass for simplicity. 
Note that the product P+P- = M x gives the invari- 
ant mass of the hadronic state. Hence the charmed 
background is eliminated as long as the phase space 
P+P- > M|, is cut out. This is achieved in various 
ways, e.g. by imposing restrictions on the leptonic in- 
variant mass q 2 = (Mb — P+)(Mg — P-)i or Ei, or 
P+, or Mx, or combinations of them. The theoreti- 
cal tools used for the calculation of such partial rates 
depend on whether the range of integrations are such 
that both P + and P_ are much larger than Aqcd (op- 
erator product expansion, OPE), or P + of order Aqcd 
while P_ much larger (QCD factorization theorems in 
this "shape- function region" ) . If both P + and P_ are 
small then resonances are not sufficiently smeared and 
the calculations are unreliable due to a breakdown of 
quark-hadron duality. 



It is convenient to use a set of dimensionless vari- 
ables which we define as y = (P_ — P+)/(Mb — P+) 
and e — 1 — 2Ei/ (M B — P+ ) , for which the phase-space 
is simply < e < y < 1. The dependence of the fully 
differential decay rate on the charged lepton energy - 
or on e in the dimensionless variables - can be sepa- 
rated from the dependence on hadronic variables, and 
the differential decay rate reads [l| 



d s T,. 
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2/(1 - V) ?2 + e(y - e) J 3 



where the structure functions 1i depend on both the 
hadronic variables P+ and y and on the factorization 
scales fj,h and but not on e. It is worth noting 
that the above formula is written entirely in terms of 
hadronic quantities, so that at this step no uncertainty 
associated with our inability to determine the 6-quark 
mass m,b is introduced. 

The differential decay rate is formally independent 
of the factorization scales fih and /i^. The scale de- 
pendence of the structure functions 3^ cancels against 
the evolution factor U y (^h, Hi), which sums (dou- 
ble) logarithms of the ratio of both scales and car- 
ries also dependence on y. The "natural" scaling for 
the scales, that avoids large logarithms, is such that 
/s 2 ~ (P + P_) and fih ~ (P-), where the brackets in- 
dicate an average over the considered region in phase 
space when integrating the differential decay rate. In 
most applications the phase space surviving the ex- 
perimental cut is such that the shape-function region, 
where QCD factorization theorems hold, is included 
and dominates the event rate. However, it is desirable 
to have a description of the decay rate in the entire 
phase space for e.g. studies of detector resolution ef- 
fects. In the next part of this talk we will first address 
the theoretical issues of the differential decay rate in 
the shape-function region, and discuss the non-trivial 
transition into the OPE region. The second part deals 
with relations between semileptonic decays and the 
radiative B — > X s ^ photon spectrum, which are con- 
structed in such a way to be model-independent at 
leading power. 
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2. Factorization 

In the shape- function region there are three relevant 
scales in the problem: the hard scale ~ ~ m b , 
the intermediate scale fi^ ~ (P + P_) ~ m b AQCD, and 
the non-perturbative scale Aqcd- When systemati- 
cally separating the physics effects at these scales, the 
structures 3^ factorize into three parts correspond- 
ingly, the hard function, the jet function, and the 
shape function Q- In modern language this factor- 
ization is achieved via a series of effective-field the- 
ory calculations QCD^SCET^HQET. The match- 
ing coefficients are by construction free of infra-red 
physics and capture the physics at the hard and 
jet scale jj, L|- Symbolically we can write 3^ = 
Hi(y,fJ,h)J(y,P+,^i) ® S(Hi) + . . ., where <g> denotes 
a convolution integral, and the ellipsis denote terms 
that are suppressed by Aqcd/™& as compared with 
the leading term. It was shown explicitly in 0, 0] 
that these power corrections factorize themselves in a 
similar fashion. The major difference is, however, that 
the leading-power shape function is unique, while at 
subleading order more than one soft function appear. 

Let us return to the leading-power formula. Not 
only is the shape function universal, but in fact also 
the jet function, and we may combine the hard func- 
tions with the kinematic prefactors to write 

J3p(0) r<1 IT/ , 12 

" ' " - Uy (M B -P + fH u (y,e) (2) 



dP+dyde 192tt 

x / doj ynibJ (ymi,(P+ —U!))S(u!) . 
Jo 

In the above expressions we have neglected the scale 
dependence, which is as discussed above. The su- 
perscript ^ on the left-hand side indicates that this 
formula is the leading power contribution only. The 
hard and jet functions have been calculated at one- 
loop order in 0, Q , and recently the two-loop result 
for the jet function has become available 0. The 
shape function S(u>), on the other hand, is not calcu- 
lable in perturbation theory and must be determined 
by other means. This is where the universality of the 
shape function comes into play. Considering the well- 
measured B — > X s ~/ photon spectrum near the end- 
point of maximal photon energy E 1 , one can state a 
formula similar to equation (J2J, namely 



(0) 



1 dT 



U {M B - P+) 3 
H r ml 

x / dili mt,J(mb(P+ — u)))S(u>) 
Jo 



(3) 



For later convenience we have normalized the photon 
spectrum to the total decay rate. 

Many of the ingredients of this equation have 
already been discussed. The main difference to 



the semileptonic decay is that a) the short-distance 
physics is different leading to a different hard func- 
tion, and b) the phase space is given by P_ = Mb 
(hence y = 1) and P + — M B — 2E 7 . Taking this 
into account, U still resums logarithms of the ratio 
/J-h/ 'Hi, but no longer carries a dependence 1 on y. We 
encounter the same shape function and jet function, 
with the latter evaluated with y = 1. 

The general idea for the determination of | V u b | from 
inclusive B decays is to use the B — > X s -f photon spec- 
trum to extract the shape function, and subsequently 
use it for predictions of arbitrary decay rates by inte- 
grating the differential decay rate This is typically 
achieved by assuming a reasonable functional form for 
the shape function in the domain where & ~ Aqcd 
with adjustable parameters, and fitting the parame- 
ters to the data. The current state-of-the-art parame- 
terizations involve exponential-type, gaussian-type, or 
hyperbolic-type functions with two free parameters, 
which can be linked to the calculation of the first few 
moments of the shape function. Therefore the extrac- 
tion of the leading shape function is mainly an exper- 
imental issue, while theory contributes via the calcu- 
lation of moment constraints (where a link to other 
decay processes is established, e.g. to B — > X c l~D, 
see for example 0, 0, 0] an d scheme translations in 
[rUp. and via the explicit calculation of the radiative 
tail of the shape function. 

Beyond the leading-power approximation there are 
several corrections that one should take into ac- 
count. Non-perturbative ("hadronic") corrections are 
encoded in subleading shape functions 0, , of which 
there are three new, independent ones at tree level 
called t(uj),u(u)) and v(uj). Their functional form is 
unknown and the only information available at present 
are their first few moments at tree level. Furthermore 
there are corrections proportional to arbitrary powers 
of the ratio of P+/mb, called "kinematical" correc- 
tions, which start at order 0(a s ) and come with the 
leading shape function. 

To make the extraction of the leading shape func- 
tion from the B — » X s "f photon spectrum feasible, 
a certain combination of t(w),u(w) and v{ui) is ab- 
sorbed into the redefined leading shape function. The 
only structure surviving at subleading power is then 
proportional to (A — u)S(u>), where A = Mb — nib is 
a heavy-quark parameter. It has recently been pro- 
posed to cross-check the findings by comparing with 
B —> X c l~v decay spectra, where a shape- function re- 
gion also exists [12j |. However, so far many important 
corrections - for example subleading shape-function 
contributions - have not been included. 



lr The relation is U y (^ h ,m) = U(fi h , /i i )-y _2a r(M) l ,Mi) i where 
orfft, fii) is the integrated cusp anomalous dimension. Sec 1] 
for details. 
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3. Transition into OPE region 

The region where an operator product expansion 
applies is reached when the size of the integration do- 
main over shape functions, i.e. the maximal J+, is 
much larger than Aqcd- The key observation is that 
the factorization theorems for the shape-function re- 
gion connect with the traditional OPE calculation via 
moments of the shape functions. More generally, if 
the integration domain A in an integral of the form 
(with some function /(£)) 

alio f{Lo)S{io) (4) 

Jo 

is much larger than Aqcd, then one can perform an 
operator product expansion in the quantity Aqcd/ A 
13] with A — A — A. An interpolation between the 
shape-function region and the OPE region, and there- 
fore a description for the entire phase space, can now 
be achieved by essentially brute force. Let us demon- 
strate the procedure for the partial decay rate with a 
cut on P + , i.e. 

r - (a) = f dP+ ^- (5 > 

We will first consider the case of small A, where the 
factorization theorems are valid, and identify a certain 
contribution (with label "SF" ) . Then we will ask what 
happens to this contribution when increasing A so 
that an OPE calculation is valid (labeled "OPE"). 

The leading-power "SF" contribution feeds into 
all powers in the OPE via the moments of the shape 
function, which involve the heavy-quark parameters 
rrib, jU^., . . ., but also misses some contributions at ev- 
ery level in the power counting. We will return to 
this last point shortly. The first subleading hadronic 
"SF" corrections (from subleading shape functions), 
which are known at tree level, do not contribute to 
the leading-power "OPE" piece, but again feed into 
all subleading powers below. Similarly the second 
subleading hadronic "SF" corrections give no contri- 
bution to the leading or the first subleading "OPE" 
piece. (Of course, when expanding in 1 /m& there is no 
first-order power correction in the OPE calculation.) 
Currently there is no complete categorization of sub- 
subleading shape functions at order \jm\. However, it 
is via the above observation that at least those which 
contribute to the \jm\ piece in the OPE calculation 
can be simulated [l|. 

As mentioned earlier, there are also kinematical ra- 
diative corrections that are power suppressed in the 
shape-function region simply because they are pro- 
portional to a s (p) ■ (P + /m b ) k , k > 1. Such "SF" con- 
tributions come with the leading shape function and 
are promoted to leading power "OPE" terms when 
P+ becomes of order nib- The exact 0(a s ) kinematic 



corrections are known by comparing the factorized ex- 
pressions with the fixed-order partonic calculation in 
full QCD [3. 

By including all of the above contribution in the 
shape- function region we correctly reproduce the OPE 
result when allowing A to become large, up to order 
l/?n^ corrections; thereby achieving a reliable inter- 
polation between the two phase-space regions. 

4. Theoretical uncertainties 

Predictions for partial semileptonic decay rates can 
be broken up into the individual contributions from 
different powers according to 

r n = rW + (r£ n «+r£ adr W) 

+ (r ki„(2) +r had r (2) ) + _ ( 

where the superscript indicates both the nature of the 
correction and their order in power counting in the 
shape-function region. For the kinematical corrections 
pkm(n) k e sum f a n terms is known and can be used 
without truncating the series. We discuss the uncer- 
tainty estimates term by term in the above formula. 

The leading-power contribution has been factorized, 
as stated in (J2J, into hard, jet, and shape function. As 
mentioned earlier, the leading shape-function uncer- 
tainty is coming from experimental limitations, and is 
not accounted for in the list of theoretical uncertain- 
ties. The hard and jet function have perturbative ex- 
pansions in the strong coupling constant a s , evaluated 
at the scale jih and jii, respectively. While is for- 
mally independent of the choice of these scales, a resid- 
ual scale dependence is introduced by truncating the 
perturbative series. A variation of the hard scale \ih 
is used to estimate the unknown 2-loop contribution 
to the hard function. The uncertainty in the jet func- 
tion has been estimated by assigning ±[a s (/Xi)/7r] 2 as 
a relative error. Here no scale variation has been used 
since a change in fj,j also changes the shape function 
S(u),/j,i), which is assumed to be extracted at a fixed 
scale, which is chosen to be /Lt« = 1.5 GeV. However, 
the uncertainty from the intermediate scale is some- 
what outdated today, ever since the appearance of the 
recent 2-loop calculation of the jet function 0- 

The kinematical corrections T^ 1 " start at order a s 
and are convoluted with the leading shape function. 
They can thus be viewed as the product of subleading 
hard and jet functions. However, for them the factor- 
ized expressions have not been worked out, and the 
effects from the scales \ih and [ii have not been disen- 
tangled. Instead, the expressions are evaluated at one 
common scale /2 which is independent of, but around 
the intermediate scale. Higher-order effects are again 
estimated by a variation of /2. 

So far we have accounted for three different pertur- 
bative uncertainties resulting from the scales fj,h, (J-i, M> 
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which can be added in quadrature. Next, we turn the 
discussion toward the hadronic corrections. It is nec- 
essary to obtain an idea about the severeness of our 
ignorance of the functional form of subleading shape 
functions. To make even a central- value prediction 
we need to model t((2>), u(&>), v{<2>) while respecting the 
moment constraints from OPE calculations. The idea 
for an estimator for the associated uncertainty is to 
find many different realistic models and study the re- 
sulting deviations from the default value. This can 
be achieved by altering the functional forms of each 
of the subleading shape functions via additional func- 
tions that have vanishing first few moments. Using 
four such functions and either adding or subtracting 
them, or leaving the subleading shape function un- 
changed, already gives a set of (2-4 + l) 3 = 729 differ- 
ent models for the set of subleading shape functions. 
To be on the conservative side, we use the maximal 
deviation from the central value as an estimator for 
the subleading shape-function uncertainty. 

At order (1/mf,) 2 we expect to find many new shape 
functions. As described earlier, we only include those 
contributions that feed into the (1/mb) 2 terms of the 
OPE result and neglect even smaller contributions. 
At this level we need not worry about their precise 
form and may model them using (/x 2 /m 2 ) S(uj) and 
(X2 / ml) S (uj) . The differences between these expres- 
sions and their true functional forms can be absorbed 
into the subleading shape-function uncertainty. 

Lastly there is one non-negligible error estimate at 
third and higher order in power counting, which is the 
weak annihilation effect. This error must be included 
whenever the experimental cut includes the phase- 
space region near the origin, where P+ ~ P_ ~ Aqcd- 
A recent study has put a limit of ±1.8% on the total 
rate by analyzing CLEO data 0]. A second possi- 
bility is to cut out this region in phase space, at the 
moderate cost (few events are located in that region) 
of a smaller efficiency. We have found that this might 
improve the overall error estimate slightly. 

In summary, we may split the theoretical uncer- 
tainty into three categories: perturbative, hadronic, 
and weak annihilation. The sizes of the individual er- 
rors in each of the categories depend on the specific 
cut employed in the measurement. For example, cut- 
ting on large charged-lepton energy is typically low in 
efficiency and is significantly sensitive to subleading 
shape functions. In this example the weak annihila- 
tion and hadronic errors dominate over the perturba- 
tive one unless the cut is relaxed below 2.1 GeV. (The 
charm background starts at 2.3 GeV.) On the other 
hand, for more efficient cuts, like a cut on the hadronic 
invariant mass or on P + for example, the perturba- 
tive uncertainty is typically dominant. At present, 
the combined theoretical error on \V u b\ - excluding 
the uncertainty from the leading shape function - is 
in the neighborhood of 5%, but may be larger for cer- 
tain cuts. 



5. Comment on "Dressed Gluon 
Exponentiation" 

A different approach in calculating inclusive B- 
decay spectra recently surfaced, named Dressed Gluon 
Exponentiation (DGE) by Andersen and Gardi 0, 
17] . This computation is based on perturbation theory 
using an on-shell 6-quark state "dressed" with gluons 
instead of the hadronic B-meson state. An interesting 
consequence of this state is that the kinematic range 
extends beyond the phase space of a single on-shell 
6-quark state. This observation motivates to use the 
calculation in comparing it with real data, so as to 
judge how realistic a dressed heavy quark can mimic 
a heavy meson. 

The procedure of dressing the 6-quark with gluons 
requires a complete knowledge of the anomalous di- 
mensions of the jet and soft functions (i.e. to all orders 
in a s ), essentially because of renormalization-group 
running deep into the non-perturbative regime. Some 
limited information is available via the perturbative 
calculation to first few orders, but the remaining as- 
pects must be modeled. Such models are guided by 
the large-/?o approximation and certain renormalon 
cancellations. 

For these reasons the DGE calculation should be 
viewed as a model, and not as a rigorous QCD pre- 
diction, the latter being a systematic expansion and 
model- independent. While it might be interesting to 
test this model against experimental data, it would be 
dangerous to use it for the extraction of \V u b\, since 
the uncertainty introduced by the underlying assump- 
tions are not under control. The statement that no 
non-perturbative function is needed and that the "pre- 
diction for the spectrum depends only on a s and on 
the quark short-distance mass" [T^ | (emphasis as in 
that reference) is not quite accurate; the same could 
be said about a simple one-parameter model for the 
shape function, where the model dependence is un- 
controlled. 

Lastly it should be noted that the DGE calculation 
has not been performed to a comprehensive level com- 
parable to the one described in the previous sections, 
i.e. many power corrections that we addressed above 
have yet to find their way into the DGE framework. 

6. Model-independent relations 

The universality of the leading shape function al- 
lows for infra-red safe relations between different in- 
clusive B-decay spectra. In fact, the program outlined 
so far - using the experimental data on the B — > X s ~f 
photon spectrum to extract the leading shape func- 
tion, and subsequently plugging this function into the 
formula for the semileptonic differential decay rate - 
can be viewed as a "manually shape-function free re- 
lation". However, a more direct relation is desirable 
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because the extraction of the shape function is some- 
what cumbersome. Such relations would eliminate the 
issues arising from the parameterizing the shape func- 
tion, and different fitting procedures that avoid sen- 
sitivity to resonances. The idea is to "invert" the jet 
function acting on the shape function in (0 so that 
instead of convolving with the shape function (which 
is not calculable) one convolves with the photon spec- 
trum (which is measured) . The desired relation reads 



r„ 



WnbV I dP+ W(AP + 



1 dT s 

r„ dP + 



\Vub\ Trhc 



(6) 



where the second term r r h c collects residual power 
corrections, e.g. terms that arise from the fact that 
different combinations of subleading shape functions 
appear in the semileptonic and radiative decay rates. 
The left-hand side of the equation denotes the par- 
tial semileptonic decay rate as obtained after impos- 
ing a cut. For the determination of \V u b\ the partial 
B — > X u l~v decay rate, as well as the photon energy 
spectrum, enter as experimental inputs. The weight 
function W(AP + ) and the residual terms r r h c , on the 
other hand, are theoretical quantities. 

The main obstacle in finding an infra-red safe 
weight function is that the jet function J(p 2 ) in Q is a 
distribution, which does not "invert" easily. Although 
the jet function is universal, it cannot be eliminated 
entirely cither, because it is called with different ar- 
guments in J2J and in ®. The difference is that in 
radiative decays y = 1, while this is generally not true 
for semileptonic decays. This problem is solved by 
defining a "jet kernel" Y(fc,lni/) which is used to ex- 
tract the y-dependence from the jet function by means 
of the following equation [Tsj 



ym b n 



m(,(fi-fc) 



dp 2 J{p 2 ) = / dkY(k,\ny) 



dp 2 J(p 2 ) 



This equation defines the jet kernel to all orders 
in perturbation theory. The explicit expressions for 
Y(k, lny) have been calculated to complete 2-loop or- 
der. We now proceed in constructing the weight func- 
tion. 

The experimental cut on the B —> X u l~v phase 
space can be encoded as 



< P+ < A , 
< y<y max 

< e < e max (P+,y) 



(8) 



As written in © the weight function depends on 
A. More generally, the weight function (and also 
the power correction r rnc ) changes as the specifics of 
the cut change, that is, it depends also on the func- 
tions ?/max(-P+) and e max (P+,y) that encode the cut. 



When integrating the leading-power differential decay 
rate formula over this phase space one can show 
via simple interchanges of integrations that the weight 
function W^°\A,P + ) itself factorizes symbolically as 

(A, P+) ~ H r H u (y, e) <8> Y(k, lny) , (9) 

where this time <g> is a three-folded convolution in the 
variables e, y and k. Because the specifics of the cut, 
A, y max (P+), and e max (P+,?/) enter only as integra- 
tion limits in this structure, weight functions can be 
computed in an automated fashion. 

Before we turn the discussion to the error analy- 
sis when using relation 10 for the extraction of \V u b\, 
we have to note a feature that contrasts direct partial 
decay-rate prediction described in Section [21 There, 
the fe-quark mass enters only indirectly into the factor- 
ization formula J5J through the shape function. When 
integrating the differential decay rate to obtain a par- 
tial rate, the nib dependence is generated dynamically. 
For example, the total decay rate is proportional to 



/ dP+ (M B 



P+ \5 



S(P+) 



at tree level. The 



analogous statement is true for the B — > X s j photon 
spectrum. This is the reason for the mjj factor in the 
denominator in @ for the normalized photon spec- 
trum. Therefore the weight function W(A,P+) will 
pick up this explicit m& dependence. Previous works 
(19 . l20j used the absolute photon spectrum, which 
introduces other, more severe uncertainties. Firstly, 
without the normalization one can only extract the 
ratio |V^b| 2 /|V^*Vtb| 2 . Secondly, the radiative correc- 
tions are then unacceptably large due to large oper- 
ator mixing effects [2l(. Thirdly, event fractions in 
B —> X s j can be calculated with higher precision than 
the absolute branching ratio |l3[ . 

As an example let us consider the theoretical un- 
certainties on the partial decay rate when using re- 
lation © for a pure cut on P + . To this end we may 
pretend for the moment that the experimental data on 
the shape of the photon spectrum had no uncertainty. 
The analysis was carried out in |2lj . and reads 



T U {P+ < 0.65 GeV) 

(46.5 ±1.4 [pert] ±1.8 [hack] ± 1.8 [m b ] 
±0.8 [pars] ±2.8 [norm]) iKbPpS^ 1 



(10) 



Again, the perturbative uncertainty stems from a vari- 
ation of the renormalization scales, and the hadronic 
uncertainty from scanning over many models of sub- 
leading shape functions. While the latter is compara- 
ble to the error analysis of the direct next-to-leading 
order partial rate prediction, the perturbative error is 
significantly reduced due to the fact that the jet ker- 
nel has been evaluated to complete 2-loop order. (In 
the direct prediction we find ±2.5 [pert] |T4b| 2 ps _1 
with the 1-loop jet function.) Varying the input pa- 
rameters needed for the weight function, i.e. quark 
masses and HQET parameters, leads to the next two 
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stated errors. Finally we have to consider the norm 
of the relative photon spectrum, which impacts |U„b| 
directly. Because the shape of the photon spectrum is 
not measured over the full kinematic region but only 
above some photon energy Eq, we will need the theo- 
retical prediction for the event fraction that falls into 
that window [lflj . 
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7. Conclusions 
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The CKM matrix element \V u b\ is a fundamental 
parameter of the Standard Model, which is already 
the first motivation in pursuing its precise determina- 
tion. It also allows for testing the CKM picture of 
CV violation and for indirect New Physics searches. 
A determination with overall uncertainty around and 
below the 10% level is both experimentally and theo- 
retically challenging, but feasible. In this talk we have 
touched upon many of the theoretical issues and pos- 
sibilities, from direct predictions of partial rates us- 
ing QCD-factorization theorems to relations between 
different inclusive i?-decay spectra, which are model- 
independent at leading power. 

The theoretical error can be reduced in several 
ways. Higher-order computations allow for a reduc- 
tion of perturbative uncertainties. The recently pub- 
lished 2-loop result for the jet function yj will have 
an impact on both the direct theoretical prediction 
of partial rates - via reduced scale dependence - and 
the weight-function relations. While the weight func- 
tion itself was already known to 2-loop order at the 
intermediate scale, the improved jet function will im- 
pact the theoretical error on the event fraction in 
B — > X s "f, i.e. the norm of the photon spectrum. 
Other possibilities for improving the theoretical errors 
include the computation of higher-order hard func- 
tions; however, a complete next-to-next-to leading 
order description in renormalization-group improved 
perturbation theory would require 3-loop anomalous 
dimensions and the 4-loop cusp anomalous dimension, 
which seems currently out of reach. On the non- 
perturbative side, further studies of subleading shape 
function contributions, and improved limits on weak 
annihilation can have a similar overall impact. 

A precise determination of \V u b\ requires a variety 
of different methods and measurements. The study of 
inclusive B decays remains a very active field of re- 
search with further improvements in both theory and 
experiment. 
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